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The Weyl curvature includes the Newtonian field and an additional field, the so-called anti- 
Newtonian. In this paper, we use the Bianchi and Ricci identities to provide a set of constraints 
and propagations for the Weyl fields. The temporal evolutions of propagations manifest explicit 
solutions of gravitational waves. We see models with purely Newtonian field are inconsistent with 
relativistic models and obstruct sounding solutions. Therefore, both fields are necessary for the 
nonlocal nature and radiative solutions of gravitation. 

PACS numbers: 98.80.-k, 98.80.Jk, 47.75.-|-f 



I. INTRODUCTION 

In the theory of general relativity, one can split the Riemann curvature tensor into the Ricci tensor defined by the 
Einstein equation and the Weyl curvature tensor[l, 2, 3, 4]. Additionally, one can split the Weyl tensor into the 
electric part and the magnetic part, the so-called gravitoelectric/-niagnetic fields[")], being due to some similarity to 
electrodynamical counterparts[2, 6, 7, 8, 9]. We describe the gravitoelectric field as the tidal (Newtonian) force[!), ID], 
but the gravitomagnetic field has no Newtonian analogy, called anti-Newtonian. Nonlocal characteristics arising from 
the Weyl curvature provides a description of the Newtonian force, although the Einstein equation describes a local 
dynamics of spacctime[9, 11]. The Weyl curvature also includes an additional force: the gravitomagnetic field that is 
produced by the mass currents analogously an electric current generating a magnetic field [2]. In fact, the theory of 
general relativity predicts two main concepts: gravitomagnetic fields and gravitational waves. Gravitation similar to 
electromagnetism propagates at identical speed, that provides a sounding analysis and a radiative description of force. 
We notice the Weyl tensor encoding the tidal force, a new force by its magnetic part, and a treatment of gravitational 
waves. 

Determination of gravitational waves and gravitomagnetism (new force) is experimental tests of general 
relativity[ I 2]. Gravitational radiation of a binary system of compact objects has been proposed to be detected 
by a resonant bar[1.3] or a laser interferometer in space[14, 15], such as the LIG0[1G] and VIRG0[17]. A non-rotating 
compact object produces the standard Schwarzschild field, whereas a rotating body also generates the gravitomagnetic 
field. It has been suggested as a mechanism for the jet formation in quasars and galactic nuclei[18, 19]. The resulting 
action of the gravitomagnetic fields and of the viscous forces implies that the formation of the accretion disk into the 
equatorial plane of the central body while the jets are ejected along angular momentum vector perpendicularly to 
the equatorial plane[2, 18]. The gravitomagnetic field implies that a rotating body e.g. the Earth affects the motion 
of orbiting satellites. This effect has been recently measured using the LAGEOS I and LAGEOS II satellites [_:;.]. 
However, we may need counting some possible errors in the LAGEOS data[21]. Using two recent orbiting geodesy 
satellites (CHAMP and GRACE), it has been reported confirmation of general relativity with a total error between 
5% and 10% [22, 23, 24]. 

In this paper, we describe kinematic and dynamic equations of the Weyl curvature variables, i. e., the gravitoelectric 
field as the relativistic generalization of the tidal forces and the gravitomagnetic field in a cosmological model contain- 
ing the relativistic fluid description of matter. We use the convention based on SttG = 1 = c. We denote the round 
brackets enclosing indices for symmetrization, and the square brackets for antisymmetrization. The organization of 
this paper is as follows. In section II, we introduce the 3-1-1 covariant formalism, kinematic quantities, and dynamic 
quantities in a hydrodynamic description of matter. In § III, we obtain constraint and propagation equations for the 
Weyl fields from the Bianchi and Ricci identities. In § IV, rotation and distortion arc characterized as wave solutions. 
In Section V, we study a Newtonian model as purely gravitoelectric in an irrotational static spacctime and a perfect- 
fluid model, and an anti-Newtonian model as purely gravitomagnetic in a shearless static and perfect-fluid model. 
We see both models are generally inconsistent with relativistic models, allowing no possibility for wave solutions. 
Section VI provides a conclusion. 
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II. COVARIANT FORMALISM 



According to the pattern of classical hydrodynamics, we decompose the spacetime metric into the spatial metric and 
the instantaneous rest-space of a comoving observer. The formalism, known as the 3 + 1 covariant approach to general 
relativity[25, 26, 27, 28, 29, 30], has been used for numerous applications [10, 31, 32, 33, 34, 35]. In this approach, we 
rewrite equations governing relativistic fluid dynamics by using projected vectors and projected symmetric traceless 
tensors instead of metrics [10, 34]. 

We take a 4- velocity vector field in a given 3+1 dimensional spacetime to be a unit vector field w°Ua = — 1- 
We define a spatial metric (or projector tensor) hab = gab + UaUfc, where gab is the spacetime metric. It decomposes 
the spacetime metric into the spatial metric and the instantaneous rest-space of an observer moving with 4-velocity 
[2, 33, 3G]. We get some properties for the spatial metric 

habu'' = 0, ha^hcb = hab, ha"" = 3. (1) 

We also define the spatial alternating tensor as 

£abc = Vabcd.u'^, (2) 

where rjabcd is the spacetime alternating tensor, 

Vabcd = -Al^\S°[J\S\S^a], Sa' = gacg'\ H = ^et gab- (3) 
The covariant spacetime derivative Va is split into a covariant temporal derivative 

Ta... =u''VbT,..., (4) 

and a covariant spatial derivative 

DbT„... =/iAa"---Vdr,.... (5) 
The projected vectors and the projected symmetric traceless parts of rank-2 tensors are defined by 

^{a> = hJ'Vb, Sl^ab) = {^(a'^^h)'' " \h'^'^hab] Scd- (6) 

The equations governing these quantities involve a vector product and its generalization to rank-2 tensors: 

[V,W]a=eabcV''W^, [S,Q]a=eabcS''dQ''', (7) 
[V, S]ab = ScdiaSbfy"^, [V,S](^ab)=£cd{aSb)''V'^. (8) 

We define divergences and rotations as 

div(l/) = D-^K, [div S) a =B>' Sab, (9) 

(Curly)„ = eafccDVS icnrlS)ab = Scdia'D'Sb)" , (curl^)(a6) = Scdia^^Sb)"". (10) 

We know that Bchab = Dd^abc, hab = '^ui^a^b), and iabc = 3u[aebc]du'' , then u°-hab = -ut and u^iabc = -u^Sabc- 
From these points can also define the relativistically temporal rotations as 

[u,V]a^ -U^iabcV'', [u,S]ab = -u''£^cl(aSb)'^, [u, S](^ab) = -^^^cd^aS b)'^ ■ (H) 

The covariant spatial distortions are 

^(aVb) = ^(aVb) - \[AwV)hab, (12) 
'^{aSbc) = ^(aSbc) - |/i(a6(divS')c). (13) 

We decompose the covariant derivatives of scalars, vectors, and rank-2 tensors into irreducible components 

Va/= -/Ua+D,/, (14) 
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where 



VfcVa = - {V{a)Ub + UaUbU.y - ^QUaVb - UaabcV" - Ua[uJ, V]b^ + DaVf,, (15) 

VcS'afc = - (^S(^ab}Uc + 2U(^aSb)dU'''Uc — ^Qu[aSb)c — '^^■{aSb)'^ ^dc 

-2e,de"(a5b)V) +D,^,„ (16) 

T^aVb = \DcV'hab - \eabcCMr\V' + Vb) , (17) 



Da'S'bc — ^^'^Sd{ahb)c - |££ic(aCUliS'b)'' + B (^^ S be) ■ (18) 

We also introduce the kinematic quantities encoding the relative motion of fluids: 

VbUa = DbUa ~ UaUb, (19) 
DfcWa jQ^afc + CTafc + i^ab, (20) 

where -ita = u'^'VbUa is the relativistic acceleration vector, in the frames of instantaneously comoving observers iia = 
U{a)j O = D°Uq the rate of expansion of fluids, (Tab = 0(0^^6) — D(a"6) ~ j^afcDcw'^ a traceless symmetric tensor 
(cab = (Ta° = 0); the shcar tensor describing the rate of distortion of fluids, and ujab = Dfa^ti] a- skew-symmetric 

tensor {uJab = ^lab]^ ^a"" = 0); the vorticity tensor describing the rotation of fluids[27, 33, 37]. 
The vorticity vector Ua [-^S, 39] is defined by 

= -^Sabci^''", (21) 

where imposes ujau"" = 0, ujab^^ — and the magnitude ui'^ = ^ujab^"''' ^ 0- Accordingly, we obtain 

Wa = -^ea6cDV. (22) 



The sign convention is such that in the Newtonian theory uj = —^V x u. 

We denote the covariant shear and vorticity products of the symmetric traceless tensors as 

[(T, S]a = SabcO'^'dS'^'^, [llI , S] ^ab) — £cd{a S b)'''^'^ ■ (23) 

The energy density and pressure of fluids are encoded in the dynamic quantities, which generally have the contri- 
butions from the energy flux and anisotropic pressure: 

Tab = pUaUb + phab + '2.q(aUb) + T^ab, (24) 
gaU" = 0, Tr"-a=0, T^ab=T^(ab)i TTafclt'' = 0, (25) 

where p = TabU°-u^ is the relativistic energy density relative to m", p = ^Tabh"'^ the pressure, Qa = —T(^a)bU^ = 
—ha'^Tcbu'' the energy flux relative to u", and nab = ?^{ab> = Tcdh"^ {au'^ b) = {h'^{au'^b) ~ ^habh'^'^) Ted the traceless 
anisotropic stress. Imposing = nab = 0, we get the solution of a perfect fluid with Tab = pUaUb +phab- In addition 
p = gives the pressure-free matter or dust solution[27, 33, 37]. 

III. COSMOLOGICAL FIELD EQUATIONS 

In the theory of general relativity, we describe the local nature of gravitational field nearby matter as an algebraic 
relation between the Ricci curvature and the matter fields, i.e., the Einstein field equations: 

Rab = Tab - \Tgab, (26) 

where Rab is the Ricci curvature. Tab the energy-momentum of the matter fields, and T = Tc^ the trace of the 
energy-momentum tensor. 
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The successive contractions of Eq. (26) on usage of Eq. (24) lead to a set of relations: 

RabWy!" = \{p + 3p), hjRbcW = -Qa, ha'^h'^Rcd ^\{P- P)Kb + TTah, (27) 

R = Ra'', T ^Ta"" ^ -p + 3p, R=-T, (28) 
where R is the Ricci scalar. The Ricci curvature is derived from the once contracted Riemann curvature tensor: 

Rab — R acb- 

The Riemann tensor is split into symmetric (massless) traceless C'abcd and traceful massive Mabcd parts: 

Rabcd ~ Cabcd + Mabcd- (29) 

The symmetric traceless part of the Riemann curvature is called the Wcyl conformal curvature with the following 
properties: 

Cabcd = C[ab][cd]i C'^ bca = ^ = C a[bcd]- (30) 

The nonlocal (long-range) fields, the parts of the curvature not directly determined locally by matter, are given by the 
Weyl curvature; propagating the Newtonian (and anti-Newtonian) forces and gravitational waves. It can be shown 
that the Weyl tensor Cabcd is irreducibly split into the Newtonian C^^^^^ and the anti-Newtonian C^f,^^ parts: 

Cabcd = C^6cd + Cabldi (31) 

CS'cd = 4{u[>,-f (32) 

C„^N^ = 2eabeUicHdf + 2ecdcU[aHb]\ (33) 

where Eab = Cacbdu'^u'^ is the gravitoelectric field and Hab = ^£acdC'^'^beU^ the gravitomagnetic field. They are 
spacelike and traceless symmetric. 

The traceful massive part of the Riemann curvature consists of the matter fields and the characteristics of local 
interactions with matter 

M'^'cd = l{p + 3pW''u[,h''^d] + Iph^ch'd] 

- 2u^'^h'-\,qd] - 2u[,/i["rf]g''l - 2u["u[,7r''lrf] + 2/i["[,7r''lrf], (34) 

Therefore, the Weyl curvature is linked to the matter fields through the Riemann curvature. 

A. Dynamic Formulas 

To provide equations governing relativistic dynamics of matter, we use the Bianchi identities 

V[ei?ah]cd - 0. (35) 

On substituting Eq. (29) into Eq. (35), we get the dynamic formula for the Weyl conformal curvature[3, 40, 4f]: 

VCabcd = --'^[aiRb]c - \gb\cR) = ~^[a{Tb\c - \gb\cTd''') = Jabc- (36) 

On decomposing Eq. (36) along and orthogonal to a 4- velocity vector, we obtain constraint {C^'^a) and propagation 
(P^^^ah) equations of the Weyl fields in a form analogous to the Maxwell equations[f 0, 42, 43, 44]: 

C\ = {diyE)a - SLU^'Hab - [a, H]a - \DaP + \Qqa 

5CTa6g'' + |[c^,g]a + 5(div7r)a = 0, (37) 



C^a = {6:iYH)a+iuj''Eab+[<7,E]a+U0a{p+p) 

+ \cnY\{q)a + i[cr, 7r]a - WlTab = 0, (38) 
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P^ab = CUT\{H)ab + 2[u,H](^ab) ^ E(ab) - QEab + [u} , E] ^ab) 
+ 3<Tc{aEb)'' - ^(7ab{p+P) ^ ^^{aQb} ~ U^aQb) 

-^'^iab) - ^©Tl'afc + 5[^,7r](ah) - ^<^''{aT^b)e = 0, (39) 

P^ab = CUrl(£;)„b + 2[U, E]^ab) + H^ab) + QHab - [t^, ] (ah) 

-3f7c(Q-fffc)'' - ^^{aQb) ~ ^[<^,q]{ab) - 5CUrl(7r)ah = 0. (40) 

The twice contracted Bianchi identities present the conservation of the total energy momentum tensor, namely 

V^T.b ^ \/'>{Rab - ^abR) = 0. (41) 

It is split into a timelike and a spacelike momentum constraints: 

C^=p+{p+p)e+ div(g) + 2Ua<?" + (TafcTT'^'' = 0, (42) 
C^a = (/O + p)Ua + DaP + (?{a) + ^^Qa + ab(t - , q\a + (div7r)a + U^Tlab = 0. (43) 

They provide the conservation law of energy- momentum, i. e., how matter determines the geometry, and describe the 
motion of matter. 

B. Kinematic Formulas 

To provide the equations of motion, we use the Ricci identities for the vector field Ua- 

2V[aV6]ltc = Rabcdv!^. (44) 

We substitute the vector field Ua from the kinematic quantities, using the Einstein equation, and separating out the 
orthogonally projected part into trace, symmetric traceless, and skew symmetric parts. We obtain constraints and 
propagations for the kinematic quantities as follows [42]: 

P3 = e + - div(M) - il'^Ua - (C^ab^"' - 'TabfT'^^) + i(p + 3p) = 0, (45) 

P\ EE UJ^^^^ + - (Ta'^LOb + icurl(u)„ = 0, (46) 

P^ab = Eab - D(a'iifc) ^ "(a^fc) + (ab) + Oc{aOb)'' + ab^ + ^ {a^ b) ^ ab = 0. (47) 

Eq. (45), called the Raychaudhuri propagation formula, is the basic equation of gravitational attraction [45]. In Eq. 
(46), the evolution of vorticity is conserved by the rotation of acceleration. Eq. (47) shows that the gravitoelectric 
field is propagated in shear, vorticity, acceleration, and anisotropic stress. 
The Ricci identities also provide a set of constraints: 

= div(w) - LUaii" = 0, (48) 
C\ = |D„e - (div(T)a + curl(t^)„ + 2[u, Lo]a -Qa^O, (49) 

C^at = Hab - CUrl(o-)at + D^a^fc) + 2u(at^fc) = 0. (50) 

Eq. (48) presents the divergence of vorticity. Eq. (49) links the divergence of shear to the rotation of vorticity. Eq. 
(50) characterizes the gravitomagnetic field as the distortion of vorticity and the rotation of shear. 
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IV. GRAVITATIONAL WAVES 

We now obtain the temporal evolution of the dynamic propagations (P^'^ab) in a perfect-fluid model (g" — Tfat = 0): 

C\ = (div^;), - SLO^'Hab - [a, H]a - iDaP = 0, (51) 

C\ = {dwH)a + Sio'Eab + K E]a + iUa{p + p) = 0, (52) 

P^ab = CUl-\{H)ab + 2[u,H](^ab} - E/^ab) - ^Eab + [uJ,E]^ab} 

+3a^^aEbf - ^(Jabip + p) (53) 

P\b = Cm\{E)ab + 2[u,E](^ab) + Hi^ab) + QHab - [w, i?] (ab> - 3(T,(,Ffa)= = 0. (54) 



To the first order, the evolution of propagation is 

P\b = ^''Eab - El^ab) - iD(aC\) - feP^fc + CUrl(p2)^, 
— ^Q^Eab — T^OE{ab} — ^Eab ~ QE^I^a^ b)"^ ~ (^cd.E^''' (Tab 

1 TJ^cd ^ _ ^cd _ 1^ TP C d , ^ C ■ d 

+-tj <Tca<Tbd—Cr CFc(ai^b)d + £cd(ai^b) ^ + £cd(ai^b) ^ 

+|e [w, £;] (ab) + 4ea,(, + e.^iaEl^u:'^ + 3(7,(, P,,) = 

+3o-c(a£^b>' - 2curl([M,£;])ab - iD(at^"i?b>c 

-|D(„ [(T,i?]fc) + fe[u,i7](„6) +curl([cj,i/])ab 
+3curl (a,(ai?b>') - fTe"ecd(aD"i?6)'^ + 2e,d(aHb{'u'^ 
+ 2ecd(aHbfu'^ + 2ecd(aHbfu'' - ^(^ab{p + p) 

-\Qaab{p + p) = 0. (55) 

We neglect products of kinematic quantities with respect to the undisturbed metrics (unexpansive static spacetime). 
We can also prevent the perturbations that are merely associated with coordinate transformation, since they have no 
physical significance. In free space, we get 

P\b = ^""Eab - Ei^ab) - |D(a b) - leP\b + cnrl{P'')ab = 0. (56) 

To be consistent with Eqs. (51)~(54), 'D^Eab — E(^ab) has to vanish. Similarly, the evolution of P^ab shows that 
D^Hab — ^f(ab) = 0. The evolutions reflect that the divergence-less and non-vanishing rotation of the Weyl fields are 
necessary conditions for gravitational waves: 

(divP)„ =(divif), = 0, curl(i;),b ^0^ cui\{H)ab. (57) 

Indeed, the rotation of the Weyl fields characterizes the wave solutions. The gravitomagnetic field is explicitly 
important to describe the gravitational waves, and is comparable with the Maxwell fields. 
We use Eq. (18) to provide two more constraints: 

C\bc = ^a Ebc - J^iaEbc) - lr)''Ediahb)c + |edc(aCUrl(S)b) = 0, (58) 
C\bc = BaHbc - ^i^aHbc) - |D'^i/rf(, hb)c + f ede(aCUrl(i/)fc)'^ = 0. (59) 

To the first order, divergence of Eq. (58) is 

^''C\bc = B^Ebc - B^'B^a Ebc) - f D^D'^Prf^, h 

+ \B''edcaCnv\{E)b'' + \ed,b^''cnr\{E)/ = 0. (60) 
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On substituting Eq. (54), it becomes 

-9|D°e''c(a^^b)d + |D''£'*c(a[t^,-ff]{h)d> 

+2DV,(aa,(,)Hd)^ = 0. (61) 
We abandon products of kinematic quantities in the undisturbed metrics: 

+ p''e^,^aP\)d - |curl(H),b = 0. (62) 
To linearized order, we get curl(S'a;))' = curl^at. Using the later point and the evolution of Eq. (53), we obtain: 



-|D(,C7\) + §D'^e%p2^)rf - lP\k = 0. (63) 

The result can be compared to the wave solution (56). Without the distortion parts, it is inconsistent with a generic 
description of wave. Distortion of the gravitoelectric field (D^(j£^f,c)) must not vanish to provide the wave solution. 
We also obtain the similar condition for the gravitomagnetic field. In free space, the divergence of the Weyl fields, 
determined by the matter, must be free. The temporal evolution decides that the rotation of the Weyl fields must be 
non-zero. Now, the distortion provides another condition to characterize the evolution of the Weyl fields: 

D(a£^6c) 7^0^D(„Ff„). (64) 

The existence of rotation and distortion is necessary condition to maintain the wave solutions. 



V. NEWTONIAN AND ANTI-NEWTONIAN FIELDS 



We can associate a Newtonian model with purely gravitoelectric (Hab = 0). Without the gravitomagnetism, the 
nonlocal nature of the Newtonian force cannot be retrieved from relativistic models. It also excludes gravitational 
waves. The Newtonian model is a limited model to show the characteristics of the gravitoelectric. We can also 
consider an anti-Newtonian model; a model with purely gravitomagnetic {Eab = 0). The anti-Newtonian model 
obstructs sounding solutions. In [43], it has been proven that the anti-Newtonian model shall include either shear or 
vorticity. 

A. Newtonian Model 

Let us consider the Newtonian model {Hab = 0) in an irrotational static spacctime {uja ~ iia ~ 0) and a perfect-fluid 
model (g° = iTab = 0). The constraints and propagations shall be 

C\ = {diYE)a - iD„p = 0, C\ = [a, E]a = 0, (65) 

P^ab^-Ei^ab) -&Eab + 3ac(aEbf - ^(Jabip + P) =0, 

(66) 

P\b = cnT\iE)ab = 0, 

= |D„e - (divd), = 0, C\b = -cnTl{a)ab^O. (67) 

To the first order, divergence and evolution of Eq. (66b) are 

B''P\b - i£a6cD^(D^S='^) + lQ[a,E]a - <yab[cy,E]'' 

= icurl(Ci)a + ieC\ - aa'Ch + ^UaP, (68) 
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P\b = -\Qcnr\{E)ab - (J,^e,d(a^^Ek)'' + curl(i;)^j, 

- \{p + p)C\b-cuT\{P^)ab + icnY\{a,(,Ek)^). (69) 

The last parameter (^Wap) in Eq- (68) vanishes because of irrotational condition. Eq. (68) then conserves the 
constraints. Eq. (69) must be consistent with Eqs. (65) and (67). Thus, the last parameters in Eq. (69) has to 
vanish: 

curl(a,(,£;b)^)=0. (70) 

It is a necessary condition for the consistent evolution of propagation. This condition is satisfied with irrotational 
product of gravitoelectric and shear, but it is a complete contrast to Eq. (65b). Thus, the Newtonian model is 
generally inconsistent with generic relativistic models. Moreover, the temporal evolution of propagation shows no 
wave solutions. 

1. Newtonian Limit 

The Newtonian model obstructs wave solution, due to the instantaneous interaction. Following [28, 29, 30], we 
consider a model whose action propagates at infinite speed (c — > c»). This is compatible with lim Eab = Eab{t)\ocn 

c — >oo 

where Eab{'t)\oo is an arbitrary function of time. 
We define the Newtonian potential as 

i;,fc = D(,Db)$ = D,Df,$-i/i,fcD2$. (71) 

On substituting into Eq. (65a), we get 

C\ = D„D2$ - iD''/i„bD2$ - iD,p = 0. (72) 

In a spatial infinity, we obtain the Poisson equation of the Newtonian potential: 

= - ip = 0. (73) 

Eq. (65a) generalizes the gravitoelectric as the Newtonian force in the gradient of the relativistic energy density. 
Moreover, Eq. (66a) gives 

P\b = -DaDfo* + \habD^^ - eDaD,,$ + \{hab + e/lah)D'$ 

+3a,(,D6)D'=$ - a,<Ji,)^D2$ - \aab{p + p) = 0. (74) 
In the Newtonian theory, we could not find the temporal evolution of the Newtonian potential. 

2. Acceleration Potential 

In an irrotational spacetime, Eq. (46) is 

P\ = \cnv\{u)a = 0. 

It introduces a scalar potential: 

ila = Da*, 

where $ is the acceleration potential. This scalar potential corresponds to the Newtonian potential 
Newtonian model, the linearized acceleration is characterized as the acceleration potential. 



(75) 
(76) 

. In the irrotational 
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B. Anti-Newtonian Model 

Let us consider the anti-Newtonian model {Eab = 0) in a shearless static spacetime (ua = Ua — 0) and a perfect-fLuid 
model ((7" = TTab = 0). The constraints and propagations shall be 

C\ = -3iU>'Hab - ^BaP = 0, C\, = (divi?), + Uaip + p) ^ 0, (77) 
P\b = cur\{H)ab = 0, P\b = H(ab) + QHab - [i^ , H] l^ab) = 0, (78) 
= fD,e + curl(Lj), = 0, C\b = Hab + 'Dia^b)+'2u^a^b)=0- (79) 

To linearized order, divergence and evolution of Eq. (78a) are 

= ^eab'D'C\ - i(p + p)C^a + Up + p)^-^' (80) 

P\b = -^ecm-l{H)ab + cm-\{H)ab 

= -|epiah + curl(p2),fc + curl([c^,i/])(,b>- (81) 

Eq. (80) is consistent only in the spacetime being free from either the gravitational mass and pressure or the gradient 
of expansion. According to Eqs. (77) and (79), the last term in Eq. (81) has to vanish: 

curl([c^,iJ])(,fc) =0. (82) 

It is a necessary condition for the consistent evolution of propagation. This condition is satisfied with irrotational 
vorticity products of gravitomagnetic, but it is not consistent with Eq. (79b): 

£^d{aC^b)c^^ - J^bC^a " jUJaC^b ~ J^b[i^,^]a - J^a[L^,Uj]b 

+ luJbBaQ + ^UJa^bQ - e''daU(b^c)^'^ " dbUia^^ c)^'^ = 0. (83) 

Thus, the anti-Newtonian model is generally inconsistent with relativistic models. Furthermore, there is not a possi- 
bility of gravitational waves. 

1, Vorticity Potential 

In an unexpansive spacetime, Eq. (79a) takes the following form: 

C\ = curl(a;), = 0. (84) 

It defines a vorticity scalar potential ^' as 

LOa - Da^-. (85) 
In the unexpansive anti-Newtonian model, the linearized vorticity is characterized as the vorticity potential. 

2. Anti- Newtonian Limit 

We may consider a gravitomagnetic model whose action propagates at infinite speed. Let us define the anti- 
Newtonian potential as 

i/ab = D(,Db)* = DaDb*-i;iafcD2*. (86) 

We substitute Eqs. (85) and (86) into Eq. (77b): 

C\ = D,d2* - iD^;i„fcD2* + D,«'(p + p) = 0. (87) 
In a spatial infinity, we derive the Helmholtz equation: 

C2 = d2* + Kp + p)* = 0. 
Eq. (77b) associates the gravitomagnetic with the angular momentum LJa{p + p)- 
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VI. CONCLUSION 



The Weyl curvature tensor describes the nonlocal long-range interactions as enabling gravitational act at a distance 
(tidal forces and gravitational waves). The gravitoelectric field is described as the relativistic generalization of the 
tidal (Newtonian) force. However, the gravitomagnetic (anti-Newtonian) force has no Newtonian analogue. We 
have no expression similar to Eab in the Newtonian theory. This difference arises from the instantaneous action 
in the Newtonian theory, which excludes a sounding solution. In § IV, the rotation and distortion of the Weyl 
fields characterize the gravitational wave. The gravitomagnetism is necessary to maintain the gravitational wave. In 
relativistic models, the Newtonian force is also inconsistent without the magnetic part of the Weyl curvature. 
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